Traffic of cytoskeletal motors with disordered attachment rates 
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Motivated by experimental results on the interplay between molecular motors and tau proteins, 
we extend lattice-based models of intracellular transport to include a second species of particle 
which locally influences the motor-filament attachment rate. We consider various exactly solvable 
limits of a stochastic multi-particle model before focusing on the low-motor-density regime. Here, an 
approximate treatment based on the random walk behaviour of single motors gives good quantitative 
agreement with simulation results for the tau-dependence of the motor current. Finally, we discuss 
the possible physiological implications of our results. 

PACS numbers: 05.40.-a, 87.16.A-, 87.16.Nn, 87.16.Wd 



I. INTRODUCTION 

Active transport on a microscopic scale is one of the 
most important features of living organisms [l[ . Trans- 
port of this kind is responsible for the functionality of 
eucaryotic cells and ultimately, via contraction of mus- 
cle cells, for the motion of organisms. Unsurprisingly, 
perturbations of the intracellular transport can lead to 
pathological conditions; for example, the role of molecu- 
lar motors in sensory defects is discussed in Q . Another 
example is Alzheimer's disease in which the breakdovifn 
of axonal transport, taking place before plaque forma- 
tion can be observed, has been reported for transgenic 
mice 0. In particular, there is experimental evidence 
that an over-expression of the linker-protein tau may re- 
duce the tubulin affinity of kincsin motors 0, Q . There- 
fore the investigation of the basic mechanisms of intra- 
cellular transport is not merely of theoretical interest. 

Studies of molecular motors discuss, e.g., the neces- 
sary prerequisites for the directed stochastic motion of 
single motors or the cooperative dynamics of many mo- 
tor proteins attached to a single cargo From a 
theoretical point of view these cases can be considered 
as single self-driven particles. Although the characteri- 
zation of the individual motions is of high relevance for 
many biological and biophysical systems, often intracel- 
lular transport is carried out by several self-driven par- 
ticles. Systems of this kind show generic many-particle 
behaviour [loj . Focusing on the many- particle features 
of intracellular transport, it is sufficient to consider an ef- 
fective motion for the molecular motors on the filaments. 
Molecular motors move generically unidirectionally and 
stepwise [ll|. The size of the steps is, in general, load 



dependent and given by multiples of the structural unit 
length of the filament. Another important feature of the 
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motor dynamics is its stochastic nature, i.e., binding and 
unbinding as well as the movements on the filament are 
random events which are described by rates. 

These properties of the dynamics suggest that intra- 
cellular transport can be described by models that are 
in close analogy to the asymmetric exclusion process 
(ASEP) — a one-dimensional lattice model for interact- 
ing self-driven particles [l^. Compared to the ASEP, 
one has to account for the finite "run length" of molec- 
ular motors, i.e., the absence of particle conservation on 
the track. Various models of this kind have been pro- 
posed, see e.g. , [3l , which show particularly interest- 
ing many-particle effects on open lattices. The absorp- 
tion and desorption of particles leads to the formation 
of localized high and low density domains. This mod- 
elling approach is very flexible and can be easily gener- 
alized, for example, in order to explain the results of in 
vitro experiments [lO[. For recent reviews treating sta- 
tistical mechanics descriptions of intracellular traffic, see 
e.g., EMI. 

In the present work we discuss the influence of a sec- 
ond type of particle which alters the attachment rates of 
the self-driven particles in different one-dimensional en- 
vironments. The model under investigation aims to re- 
produce the generic features of the kinesin dynamics on 
microtubules for different concentrations of tau. In addi- 
tion to shedding light on this speciflc biological situation 
our study contributes to the understanding of disorder in 
driven many-particle systems. 

The rest of the paper is organized as follows. In Sec.HIl 
we review the lattice-based model for intracellular trans- 
port introduced by Klumpp and Lipowsky fisj and dis- 
cuss simulation results obtained by generalizing it to in- 
clude the effect of tau. Then, to gain further insight, 
we introduce in Sec. IIIII a somewhat simplified model 
which retains the key features of the problem but is more 
amenable to analytical treatment; we also discuss some 
limits in which this model is exactly solvable. In Sec. IIVI 
we address the biologically important question of how the 
mean current and its fluctuations depend on the concen- 
tration of tau. A random walk treatment for single mo- 



tors leads to predictions for the fow-motor-density regime 
of the multi-particle model which are in good agreement 
with simulation. Finally, in Sec. |Vl we conclude by dis- 
cussing the possible physiological relevance of our results 
and potential generalizations to more realistic models. 

II. FROM BIOLOGY TO SIMULATION 

In this section we outline how the complex biological 
system of interest can be represented by a lattice-based 
model amenable to simulation. Specifically, we build on 
the work of Klumpp and Lipowsky [isj , generalizing their 
model to include tau proteins and presenting specimen 
simulation results. This is in the spirit of a now well- 
established approach to complex systems whereby one 
considers simple models in order to help understand phe- 
nomena and mechanisms which may also occur in more 
realistic situations. Indeed, we will pursue this approach 
further in the following section in which we discuss an 
even more simplified model allowing for analytical treat- 
ment. 



A. Model of Klumpp and Lipowsky 

As discussed in the introduction, motor proteins carry 
out directed walks (i.e., active transport) on the filaments 
of the cytoskeleton. In the context of neuronal transport, 
we are particularly interested in the motion of kinesin 
along axon microtubules (composed of repeating tubulin 
subunits) [l^ . The stepwise structure of such movement 
naturally suggests using a one-dimensional (ID) lattice 
model with lattice-spacing equal to the step size [ll| . We 
assume that the size of the kinesin molecules is compa- 
rable to this unit lattice spacing; the mutual exclusion of 
motors is then imposed by preventing double occupancy 
of any one site. An obvious starting point for modelling 
the dynamics is the prototypical asymmetric simple ex- 
clusion process in which particles hop to vacant nearest 
neighbours with some preferred direction (for a review of 
this exactly-solvable model, see [2fi|)- In fact "backward" 
steps for kinesin are extremely rare so it suffices to con- 
sider the totally asymmetric version in which particles 
can only hop in one direction. 

Since kinesin motors typically detach from the micro- 
tubule after a certain number of steps ( "run length" ) 
and diffuse in the surrounding medium before reattach- 
ing, one needs to extend the standard ASEP picture. A 
simple choice is to impose the same lattice structure and 
exclusion rules off the filament but allow the molecules 
to move symmetrically there. The association onto and 
dissociation from the track is governed by explicit rates 
(typically the dissociation rate is comparatively small so 
that a single motor can make many steps before detach- 
ing). 

Klumpp and Lipowsky embed this lattice model in 
a cylindrical geometry with directed motion along a 
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FIG. 1: (Color online) Graphical representation (not to scale) 
of the key ingredients in the simulation model (after [l^ but 
with the addition of tau). The cylindrical geometry can be 
considered as a simplified representation of a nerve cell where 
the track along the symmetry axis is the axon microtubule. 
Our simulations are based on a cubic lattice with periodic 
boundary conditions (not shown). Two-headed kinesin mo- 
tors (with attached cargoes) diffuse in the bulk and undergo 
directed transport along the track. The tau proteins (repre- 
sented here by small hexagons) also diffuse in the bulk but do 
not move along the track. 



ID track located at the symmetry axis of the cylinder 
and undirected motion (diffusion) elsewhere. Of course, 
this is a considerable simplification compared to a real 
nerve cell which contains a bundle of parallel micro- 
tubules. However, it is expected to reproduce qualita- 
tively the main features of unidirectional transport (and 
could readily be computationally extended to the case of 
several parallel tracks). Another important question con- 
cerns the boundary conditions at the ends of the cylinder. 
In [l^ both periodic and open boundary conditions are 
considered; in the former case, the current-density rela- 
tion is obtained and, in the latter, the phase diagram. 

B. Effect of adding tau 

In this subsection we generalize the above-described 
model of Klumpp and Lipowsky to include the effect of 
tau proteins. The interested reader is referred to [2l|, [22| 
for details of the structure of tau and its interaction with 
tubulin. Here we merely note that tau molecules "deco- 
rate" the microtubule (indeed they arc believed to play 
a role in stabilizing it) but do not undergo active trans- 
port along it. We can therefore incorporate them in our 
model by introducing a second species of particle which 
is allowed to diffuse in the bulk, and absorb/desorb with 
certain rates but not to move along the track itself. Fig- 
ure [1] summarizes the various processes now incorporated 
in the model. In the present work we restrict ourselves 
to periodic boundary conditions. Although one expects 
that real nerve cells are better described by open (or half- 
open) boundary conditions cf . [23j , studying the simpler 
case of periodic boundary conditions already gives infor- 
mation about bulk effects. We further justify our choice 
by remarking that the finite run lengths of the molecular 
motors might be expected to reduce the dependence of 
currents on the filament boundary conditions. 

To complete our description of the model we must spec- 
ify the details of the tau interactions. Firstly we note that 
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tau molecules are relatively small compared to the motor- 
cargo complex. Single motor in vitro experiments 
also indicate that the presence of tau does not affect the 
speed of kinesins along microtubules or their run length 
(probably because kinesin and tau have different binding 
sites [l^l) but does lead to a significant reduction in the 
attachment rate. 

In the present simulation model, wc therefore neglect 
tau-tau and tau- kinesin exclusion effects (i.e., each unit 
lattice step can contain one or more tau particles, regard- 
less of the presence or absence of a kinesin motor) . How- 
ever, the crucial point is the effect of tau on the kinesin 
dynamics. In vivo experiments Q are complicated by the 
fact that many motors can be attached to a single cargo 
but also suggest that one effect of tau is a decrease in the 
motor attachment rate (especially for kinesin). Hence, 
in the present model, when one or more tau molecules 
occupies a site on the track the rate for absorption of a 
kinesin is reduced.^ 

The description above furnishes a model which retains 
the key features of the biological system but can be im- 
plemented simply in numerical simulation. In principle, 
one can try to match the parameters of the model (lat- 
tice spacing, rates, etc.) to those known from experi- 
ments in order to obtain quantitatively meaningful re- 
sults. We carried out simulations using relative values 
for the kinesin rates based on data for the case without 
tau in [l^.^ Specifically an unbound kinesin could hop 
to a neighbouring empty filament site with rate 0.0083 
which was identical to the hop rate between a pair of 
empty bulk sites (i.e., the binding probability was one). 
The rate for a forward step on the filament was 0.0099 
and the desorption rate from the filament to each of the 
four adjacent unbound lattice sites was 1.7 x 10~^. The 
tau molecules were (arbitrarily) set to diffuse twice as 
fast as the motors and had absorption and desorption 
rates, 0.03 and 2.5 x 10~^ respectively, which are consis- 
tent with data from FRAP (Fluorescence Recovery After 
Photobleaching) experiments [2^ . Note the high affinity 
of tau for the microtubule which leads to significant tau 
coverage on the filament even for low global densities. 

In order to illustrate the possible effect of tau on ki- 
nesin binding and transport, we reduced the kinesin ad- 
sorption rate by a factor of one hundred in the presence 
of tau. This value arguably overestimates the impact of 
tau in real systems. (In |5| it was shown that the bind- 
ing rates depend on the ratio between tau and kinesin 



Note that some aspects of our model are similar spirit to the dis- 
cussion of crowding by another molecular species which already 
appears in [23l . However, the specifics are rather different — 
in particular, in that work, there is exclusion of motors by the 
"obstacles" and the obstacles can also be actively unbound by 
processive motors. 

The rates are equivalent to probabilities per timestep; for com- 
parison with real data, normalization of the timcscale would then 
be achieved by associating the length of each timestep with a 
physical time unit. 
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FIG. 2: (Color online) Average current (jb) along the micro- 
tubule filament (normalized by the velocity of bound motors 
Vb) as a function of kinesin density Nmo/L for different tau 
densities Qt- Note that here for consistency with [l^i the mo- 
tor density is measured in number of motors per unit lattice 
length, along the axonal cylinder — this can trivially be con- 
verted to the number per unit volume using the dimensions 
of the system (length 1000 lattice sites, radius 25). Simulation 
data is averaged over 50 histories and approximately 2 x 10*' 
time steps. Lines are provided as an aid to the eye. 



concentration; for the highest concentration of tau used 
in that experiment the binding rates were reduced by a 
factor of five.) However, we expect a stronger influence 
of tau in living cells since the typical run length of the 
molecular motors must be shorter than in in vitro as- 
says, due to the bidirectionality of the transport and the 
dynamics of the filament. 

Figure [2] shows simulation results for the current along 
the track versus bound kinesin density for different val- 
ues of global tau density. One sees from these results 
that increasing the density of tau reduces the height of 
the maximum and shifts it to higher kinesin densities. 
The average current is reduced for low densities of kinesin 
but increased for high densities. Furthermore, the sim- 
ulation results cannot be reproduced by a simple mean- 
field approximation which indicates the significance of 
correlation effects. Indeed explicit measurements of the 
kinesin pair-correlation functions show the existence of 
long-range correlations. To explore the tau-induced ef- 
fects analytically, in the next section we present and dis- 
cuss an even simpler model. 



III. MINIMAL MODEL 

In order to elucidate the mechanisms leading to the 
features observed in the simulations described above, we 
now introduce a simpler model which can be treated an- 
alytically in certain limits. Specifically, instead of the 
cylindrical geometry we consider a two-lane system (anal- 
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FIG. 3: Schematic representation of simplified model. The 
larger grey circles denote K-particles and the smaller black 
circles are r-particles. Arrows show possible moves with as- 
sociated rates. The model has periodic boundary conditions. 



Ogous to the "two-state approximation" of [l8|) with 
somewhat simplified effective dynamics for the tau par- 
ticles. 

Note that this is now a stylized "toy" model designed 
to yield increased understanding but not quantitative 
agreement with experimental results. We further remark 
that simple two-lane models have previously been used as 
crude representations of vehicular traffic and are of more 
general interest in building up our knowledge of non- 
equilibrium statistical physics, see e.g., |26l - [29| . In the 
following analysis we mainly use the language of physics 
("lattice sites", "particles", etc.) whilst keeping in mind 
the original biological motivation. The physiological sig- 
nificance of our results will then be discussed in more 
detail in Sec. El 



A. Model definition 

We consider a two-lane lattice gas model with periodic 
boundary conditions, as shown schematically in Fig. [31 
Here the lower lane represents a microtubule (MT) fila- 
ment and the upper one the surroundings. The model 
contains two species of particles representing kinesin mo- 
tors (K) and tau proteins (r). We denote the occupation 
number for K-particles on the ith site in the lower or 
upper lanes by bi or Ui respectively (where the letters 
denote "bound" and "unbound" in allusion to the bio- 
logical context). The r-particles are found only in the 
lower lane (see below) with occupation number r,;. The 
total number of each kind of particle is conserved and 
a hard-core interaction prevents more than one of each 
species on a given site, i.e, the site occupation numbers 
can take only the values or 1. The dynamics of the 
model (defined in continuous-time) are motivated by the 
biological picture, as will be explained below. 

We first make some remarks on the behaviour of the 
T-particles. After the discussion of Sec. |TT1 the exclu- 
sion condition may seem very unrealistic. However, for 
the kinesin transport the relevant question is whether 
a given location on the microtubule contains at least 
one tau molecule (in our simplified description the num- 
ber of such molecules is irrelevant). It therefore seems 



reasonable to reduce the state-space and consider only 
"tau-occupied" and "tau-empty" sites, in other words 
the presence or absence of an effective r-particle. Sim- 
ilarly, we argue that tau desorption, diffusion and ab- 
sorption can be represented by an effective undirected 
motion along the track. Hence, in our simplified model, 
the r-particles undergo a symmetric exclusion process 
(SEP) in the lower lane. Specifically, a r-particle at site 
i in the lower lane hops randomly (after an exponentially 
distributed waiting time) to one of the nearest-neighbour 
sites I — 1 or i + 1 . The rate for each of these moves is 
if the destination site is vacant and zero otherwise. It is 
clear that the average occupation number in the steady 
state pr = (Ti) is site-independent. [Here, as throughout 
the paper, we use angular brackets to denote an average 
over stochastic histories.! 



K-particles can occupy sites in both the lower lane and 
upper lane (corresponding to kinesin motors being bound 
or unbound respectively). Along the lower lane the dy- 
namics is given by a totally asymmetric exclusion pro- 
cess with rate p for rightward hops. Diffusion of motors 
in the surroundings of the filament, is represented by a 
symmetric exclusion process in the upper lane with rate 
D. The key feature of the model (reflecting the under- 
lying biology) is that the coupling of these two processes 
depends on the local density of r-particles. Specifically, 
a K-particle in the upper lane hops to a K- vacant neigh- 
bouring site in the lower lane with rate uja if the desti- 
nation site is not occupied by a r and with a (reduced) 
rate Qa if there is a r there. In contrast, a K-particle 
in the lower lane moves to a vacant neighbouring site in 
the upper lane with rate ojd, regardless of the presence 
or absence of r. 



To summarize, our model is defined by the possible 
moves (and corresponding rates) listed in Table IH^ Even 
this simple model is difficult to treat analytically since 
the K-particles experience a dynamic disorder due to the 
interaction with r. For reviews of earlier work on disorder 
in driven diffusive systems see, e.g., [33,[3l|. For a disor- 
dered asymmetric exclusion process with non-conserving 
sites, Evans et al. [s^] were able to solve the steady state 
exactly in two limits. In a similar spirit, we treat be- 
low some tractable limits of the present model as well as 
discussing more qualitatively the general case. 



^ Of course, this simplified representation would be inappropriate 
for addressing questions about the distribution and dynamics 
of the tau particles themselves. With respect to properties of 
the kinesin transport, the validity of the approach will later be 
further justified by the similarities between simulation results 
from this minimal model and from the full model of Sec. [U 

^ Note that we sometimes use the generic lo^ for statements which 
apply separately to all three inter-lane rates u)a, i^a, i^d- 
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TABLE I: Table of possible moves 



Move 



Rate 



Biological interpretation 



{n = 1, Ti±i = 0} ^ {n = 0, Ti±i = 1} 
{bi = 1, b,+i = 0} ^ {b, = 0, fe,+i = 1} 
{ui = 1, Ui±i — 0} ^ {ui = 0, Ui±i = 1} 

= l,u,^0}^ {b, = 0, = 1} 
{n = 0, 6, = 0, u, = l}^ {n = 0, = 1, u, = 0} 
{n = l,bi = 0, Ui = l} -)■ {n = 1, 6i = 1, m = 0} 



P 
D 



Effective diffusion of tau along MT 
Directed motion of kinesin along MT 
Diffusion of kinesin in surroundings 
Detachment of kinesin from MT 
Attachment of kinesin in absence of r 
Attachment of kinesin in presence of r 



B. Tractable limits 

For a fixed density of r-particles, the behaviour of the 
system is characterized by plotting the stationary current 
along the lower lane against the density of K-particles. In 
this subsection we show how this "fundamental diagram" 
can be calculated exactly in certain limits and compare 
the T-induced changes with the observations from the 
more realistic model in Sec. |TT1 



1. Pure case 

In the absence of t (i.e., pr — 0), the model is ex- 
actly solvable as shown by Klumpp and Lipowsky (isj . 
The steady state is a product state with site- independent 
K densities in lower and upper lanes ph = (bi) and 
Pu = (ui). Since this stationary state has no correla- 
tions between sites, these densities obey the mean- field 
equation 

WdP6(l - Pu) = t^aP«(l - Pb) (1) 

which corresponds to the absence of a net current be- 
tween the two lanes. Combined with the expression for 
the (known) total motor density 



this yields a quadratic equation for p^ which can be solved 
explicitly. The mean current (j) along the track is then 
obtained as 

ii) = ppbil - Pb)- (3) 

In Fig. |4]we plot {j)/p against the total motor density 
Pk] note that in the pure case the shape of this "funda- 
mental diagram" is completely determined by the ratio 
of absorption and desorption rates r = ujd/^a- 



2. Fast diffusion of r-particles 

If the r-particles diffuse infinitely fast relative to the 
movement of the K-particles then there is no possibility 
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PA- 
FIG. 4: (Color online) Current of K-particles (representing 
kinesin motors) along track versus K-density. Parameters are 
chosen so that the ratios of desorption and absorption rates in 
the presence and absence of r are given by r = u^/ua = 0.1 
and f = ujd/SJa = 0.5. Lines show (i) the exact theoretical 
result for the pure case (pr = 0) and (ii) the simple mean-field 
result for the case with pr = 0.5. Points denote simulation 
results from a system with lattice size L = 1000 averaged over 
1000 realizations and 10000 time steps. Individual cases are 
(a) pure model with parameters: Pt = 0, p = 0.2, D — 0.5, 
ujd ~ 0.05, ijJa = 0.5; (b) fast r diffusion with parameters: 
Pr = 0.5, Dr = 5.0, p = 0.2, D = 0.5, coa = 0.05, cua = 0.5, 
uja = 0.01; (c) case of small absorption/ desorption rates with 
parameters pr = 0.5, = 0.01, p = 2.0, D = 5.0, uJd = 0.05, 

U]a ~ 0.5, UJa = 0.01. 



for correlations to build up in the system and the pres- 
ence of the r-particles merely alters the average absorp- 
tion rate. Hence in the limit Dr — )■ oo (with the other 
rates all finite), one expects that the stationary state is 
still described by the simple mean- field theory (i.e., with 
uncorrelated homogeneous site densities) with an effec- 
tive absorption rate given by 

Wa, off = (1 - PT)Wa + Pr^a. (4) 

This expectation is confirmed by the simulation results 
shown in Fig. 01 Notice that the effect of the introduc- 
tion of r (with UJa < UJa) is to shift the position of the 
maximum towards higher K-densities. 
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(i) limit u)x/ D,p/ D — > oo 

(ii) limit u}x/p, D/p — > oo 

(a) large u}x/D,p/D + 

(b) large uj^Jp, D/p 




FIG. 5: (Color online) As Fig. [4] but for the quenched disor- 
der limit Dr 0. Lines show theoretical results for (i) in- 
finitely fast absorption/desorption and K-movement on track 
compared to K-diffusion in bulk and (ii) infinitely fast ab- 
sorption/desorption and K-diffusion in bulk compared to K- 
movement on track. These two cases correspond to a homoge- 
neous density profile in the lower and upper lanes respectively. 
Supporting simulation results are shown for parameter values 
(a) p-r = 0.5, D-r = O.Of, p = 2.0, D = O.Of, cj^ = 0.5, 
LOa = 5.0, Oja = O.f and (b) pr = 0.5, Dr = O.Of, p = O.Of, 
D = 5.0, ujd = 0.5, uja = 5.0, ua = o.f. 



As an aside, we note that the same mean-field solution 
is found for any value of Dt in the limit of vanishing cou- 
pling between the lanes, uix ^ Q (i.e., lUa — 0, Wa — >■ 0, 
—7- 0) with the ratios r = LOd/uJa and f = ujd/Cja con- 
stant. In this limit the time between each lane-changing 
event tends to infinity so that after each such event the K- 
density profiles in the lower and upper lane relax to the 
spatially-uncorrelated homogeneous case and we again 
recover the mean-field result with effective rate given 
by (HI (see again the simulation results in Fig. S]). 



3. Slow diffusion of r-particles 

Let us now consider the quenched disorder limit, i.e., 
Dr — > 0. In this case, the time between movement of 
T-particlcs tends to infinity and between these events 
the system relaxes to a quasi-stationary state with site- 
dependent densities of K-particles. There is a grow- 
ing body of literature discussing lattice gas models with 
quenched spatial disorder but analytical progress is diffi- 
cult even for the usual single- lane ASEP (see, e.g., [H- 
Isrj). For our two-lane model, we focus here on two lim- 
iting cases where the analysis again becomes relatively 
straightforward; comparison with simulation is shown in 
Fig. El 

Firstly, we consider taking Dr followed by 
U!x/ D,p/ D oo whilst holding the other ratios be- 
tween rates constant. Here the density profile in the 



lower lane relaxes rapidly to stationarity between each 
hopping event in the upper lane. In other words, the 
profile in the lower lane is that of a homogeneous ASEP 
(constant density pb) whereas the density in the upper 
lane is inhomogeneous and determined by the local ab- 
sorption and dcsorption rates. Specifically, there are two 
distinct unbound densities pu and pu corresponding to 
the absence and presence of r respectively. These densi- 
ties must obey the relations 



Pu) 
Pu) 



Pu{l 
Pu{l - 



- Pb) 
Pb) 



(5) 
(6) 



and the equation for total K-particle density becomes 
1 , 



Substituting ([5]) and ([6 



PK ^ 1^ 



Pb 



- (1 - Pr)Pu + PrPu] ■ 

into d?]) gives 
Pr)Pb , rprPb 



1 + (r - l)pb 1 + (f - l)pb 



(7) 



(8) 



For fixed px this is a cubic equation for pb] numerical 
solution followed by substitution into ^ yields a current 
in good agreement with simulation.^ 

Secondly we consider the case where Dr ^ and then 
uJx/p,D/p — >■ oo (again with the other ratios between 
rates held constant). Here the density in the upper lane 
relaxes to that of a homogeneous SEP (constant density 
Pu) and the density in the lower lane is inhomogeneous 
with two possible values pb and pb given by the local 
absorption and desorption rates via 



rpb{l-pu) =Pu{l-pb) 
fpb{l - p„) = pu{l - Pb)- 



These densities must satisfy 



leading to 



PK 



PK = \ [(1 - Pr)Pt + PrPb 



PbPu 



(1 - pDpu ^ 

r + {l-r)pb f+{l-r)pu 



(9) 
(10) 



(11) 



(12) 



After solving this cubic equation for pu one can obtain 
the two bound densities via dS]) and dlO]). The current 
along the lower lane is of order p and its average value 
can be calculated by using the probabilities to find a r- 
particle at either end of a given bond: 

0') = P[(l - PrfPbi^ - Pb) + Pt{1 - PT)pb{l - Pb) 

+ (1 - Pr)PrPbil - Pb) + (Pr)'pb(l - Pb)]- (13) 



Even better agreement would presumably bo obtained for Dr 
smaller relative to D but then longer simulation times would be 
required to reach the steady state. 
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FIG. 6: (Color online) Current of K-particles along track ver- 
sus total K-density per unit length for pr = 0.75, Dr ~ 0.01, 
p = 0.5, D = 1.0, LOd = 0.01, LJa = 1.0, Cja = 0.01. Crosses 
show simulation results (parameters as in Fig. 2]); lines show 
analytical results for the limiting cases of IIIIB flllll B 3l 



significant correlations were all observed in simulations 
of the more realistic model in Sec. |IT1 

From a theoretical point of view, the shape of the fun- 
damental diagram is interesting (in particular, the ef- 
fect of tau on the position and height of the maximum) . 
Physiologically however, one expects the number of ki- 
nesin motors to be approximately constant, and thus a 
more relevant question is the dependence of current on 
tau concentration at a fixed low density of motors. Of 
course, it is possible to analytically obtain this depen- 
dence in the tractable limits discussed above. For more 
"realistic" intermediate parameters, a good approximate 
description in the low K-density regime can be obtained 
by considering the random walk properties of single mo- 
tors which is the approach we choose to pursue in the 
following section. To simplify the analysis while retain- 
ing the most important features we henceforth set Dr = 
and uia = 0.^ 



IV. TAU-DEPENDENCE OF CURRENT 



Of course, on very long timescales the positions of the r- 
particles and hence the local densities of K-particles will 
change. However, the system will rapidly relax to a new 
quasi-stationary state which, on average, has the same K- 
current (as follows from conservation of r-particle num- 
ber). The prediction of (|13p is again well-supported by 
simulation even for parameter values, which for compu- 
tational efficiency, are still relatively far from the limiting 
case; see Fig. [S] 



C. Intermediate parameters 

All the limits treated in the previous subsection corre- 
spond to cases where in the (quasi-)stationary state there 
are no correlations between the occupations of different 
sites (i.e., they are situations which can be described by 
forms of mean- field theory). However, we emphasize that 
this is not expected to be the case for general choices of 
rates. In this subsection we present illustrative simula- 
tion results for intermediate parameters — see Fig. [6l 

The relative values chosen are inspired by the biological 
context although we make no attempt to quantitatively 
match our simple model to real data. Specifically, the 
bulk diffusion and absorption rate are large compared to 
the hopping rate on the track but the desorption rate is 
small. The presence of tau strongly reduces the absorp- 
tion probability (wq ^ Wa) while the effective diffusion 
rate for tau is rather small. The simulation data shown 
in Fig. ini shows that, just as in the limiting cases dis- 
cussed in IHI B 2hllH B 31 there is a shift in the position 
of the maximum towards higher K-densities (compared 
to the pure case). However, significantly, there is also a 
decrease in the height of this maximum. We remark that 
a shift in the maximum, a reduction in its height and 



A. Single-motor properties 

In this subsection we show how the properties of a sin- 
gle motor (or equivalently an ensemble of non-interacting 
motors) can be analytically understood through a ran- 
dom walk picture. This facilitates an investigation of the 
effect of tau on the average velocity and diffusion con- 
stant of the motor (or equivalently the mean and vari- 
ance of the current). Such a single-particle random walk 
approach (cf. [l3i . ,38, . ,39J ) is relevant for in vitro exper- 
iments involving individual motors but, as we shall see 
in IIVBI also gives valuable information about the low- 
motor-density limit of the interacting-motor problem. 



1. Probability distributions, large deviations 

Let us start by considering a single motor moving in 
the quenched disordered landscape defined by Fig.|3]with 
Dr = and cja = 0. We are interested in the distribution 
of the total number of steps Yi (t) made by a single motor 
in the lower lane over some time period [0,t] (i.e., the 
integrated current in time t). Now it is obvious that, if 
the motor spends a fraction x of its time in the lower 
lane, then Yi has a Poisson distribution with mean pxt, 
i.e., 

P,,v-„„.f:q?^ ,14, 



We can appeal to physical intuition and simulation results 
(cf. Fig. [FHSll to support the assertion that, with respect to the 
current of i^-particles, there is no singular behaviour at this limit. 
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or equivalently the observed time-averaged current yi 
Yi/t has the distribution 



(15) 



The full distribution of yi is then obtained by averaging 
over all possible values of x, i.e., 

Piyi,t)^ f P,{yi,t)P{x,t)dx, (16) 
Jo 

where P(x, t) is the distribution of the fraction of time 
(from a total t) spent in the lower lane. Our central task 
is thus to calculate this distribution and hence obtain the 
distribution of j/i. We note that P{yi,t) is expected to 
have a large deviation form 



P{yi,t)^e 



-te(j/i 



(17) 



and that knowledge of the large deviation function e(j/i) 
gives the asymptotic behaviour of all moments of the 
distribution. 

Before showing how to obtain P{x, t) approximately in 
the disordered case, we first digress to discuss the pure 
result corresponding to pr = 0. In this case, the upper 
and lower lanes form an effective two-state Markovian 
system with rates uja and ujd for transitions between the 
states. It can then be shown exactly that the distribution 
of the total fractional time in the lower (bound) state is 
given by 



Pix, t) = ,5(1 - X)e"(""+'^'')^"' + [LUa + LOcdP^J-x 



xP2 



where 8 is the Dirac delta function, /g and 1\ are Bessel 
functions, P° and P° are the probabilities to find the 
motor in the bound and unbound states respectively and 



X = 2{Ua + UJd)t^P^P^^x{l~x). 



Furthermore, it is trivial to show that 



pO ^ 



(19) 



(20) 



The expression ((T5)) can then be substituted into 
and numerical evaluation of the resulting integral agrees 
with the current distribution obtained by simulation (not 
shown) In the long-time limit reduces to a Gaus- 
sian [4]| . 



2. Disordered model as trapping problem 

In the disordered case (non-zero density of tau) the 
upper and lower lanes no longer form an effective two- 
state Markovian system (the average rate for escape from 



the upper lane is dependent on the occupation time, since 
for longer times the particle is more likely to be found 
trapped in long tau-decorated regions). However, for a 
two-state non-Markovian system, it can be shown [42| 
that the asymptotic distribution of occupation times is 
still Gaussian with the fraction x of time in the lower 
lane having mean 



Ph = 



and variance given by A/t where 



(21) 



(22) 



Here {{ti,)) and cr,^ (cr^) are the mean and variance of 
the sojourn time distribution in the upper (lower) lane. 
Since the rate for exiting the lower lane is always cod we 
have 



(4) - — , 



1 



(23) 



and the only remaining difficulty is to calculate the equiv- 
alent quantities for a sojourn in the upper lane. Obvi- 
ously, (tu), {t\) and hence Pb and A are functions of pr 
but for simplicity we suppress this dependence notation- 
ally. 

Now, focusing in on the upper lane, we see that a K- 
particle performs an ordinary random walk (with hop- 
ping rate D) in regions of sites occupied by tau. When 
it reaches a site without tau, it has a probability A = 
uja/{2D -\- uja) to be "trapped" and move to the lower 
lane. For several decades such one-dimensional trapping 
problems have been extensively analysed in the literature 
(see, e.g., j43l - l45| and references therein). In particular 
it can be shown (see, e.g., [1^) that the tail of the so- 
journ time distribution has a stretched exponential form 
P{tu) ^ exp(— a pr))^^^ t]/^) due to the survival of 
particles in arbitrarily large trap-free regions. However, 
we are interested in the mean and variance which are not 
much influenced by the tail of the distribution. Here we 
show how, within certain approximations, these quanti- 
ties can be calculated directly. Our approach is in the 
spirit of a continuous time- version of the method given 
in p. 

First we consider two traps separated by a "line seg- 
ment" of length as shown in Fig. [71 In other words, we 
have a tau-free site at ? = followed by — 1 sites occu- 
pied by tau, followed by another tau-free site si I = N . 
The time to reach one of the traps at either end start- 
ing from position I is a random variable Ti (N) ; it can be 
shown by standard arguments [48| that its mean is given 
by 



{Ti{N)) 



2D 



(24) 



Note that this expression is invariant under the transfor- 
mation I ^ N — I as expected from the symmetric nature 
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D 



o 



N - 1 



A : 



2D + UJa, 



FIG. 7: Schematic of a K-particle diffusing on the A'" — 1 tau- 
occupied sites between two tan- free sites (traps). This can be 
considered as a trapping reaction with trapping probabihty 

A = ijJa/{2D +ijJa). 



of the diffusion. We now make the important assumption 
that the distribution of entry points into the upper lane 
has a uniform distribution. This is expected to be a good 
approximation when the typical run length in the lower 
lane is large compared to the spacing of traps, i.e., when 



P 1 
— » . 



(25) 



Using T{N), without subscript, to denote the time to 
absorption from a random initial position we then obtain 



N-l 

1=1 

_ N'^ -1 
^ 12D 



(26) 
(27) 



Following [48|, one can also write down expressions 
(see appendix) for the mean square first passage times 
{Ti{N)^) and (T(iV)2) = (1/A^) W)- Com- 

bined with a known distribution of traps (i.e., the distri- 
bution of lengths N) , this is sufficient to obtain explicit 
expressions for the first two sojourn time moments in 
the case where the trapping probability A is unity (i.e., 
uja oo). However, for imperfect trapping one must al- 
low for the possibility that the particle is not absorbed 
on its first visit to the trap. We now discuss two possible 
approximations which allow further analytical progress: 

1 . In the "reflecting trap" or "decoupled-ring approxi- 
mation" [4^ [4^ it is assumed that a particle which 
is not absorbed by a trap is always reflected back 
into the same trap-free line segment from which it 
camc.^ This simplifies the analysis by decoupling 
the tau segments between traps but at the expense 
of introducing an uncontrolled approximation. The 
approach is clearly exact in the case of perfect traps 
(i.e., in our model for oja — > 00) but has also been 
shown to give good results for trapping probabili- 
ties smaller than unity (49| . 



2. Another approach is to place each particle which es- 
capes from a trap site on to the first site (i.e., ^ = 1) 
in a line segment with length N randomly drawn 
from the distribution of all possible lengths. Again 
this approximation is exact for uja ^ 00. Further- 
more, we might also expect it to be exact in the 
limit of zero-trapping probability (w^ 0) since in 
this case each particle explores the full distribution 
of tau line segments. 

Whereas the first approximation treats the line segments 
as unconnected, the second assumes they are all inter- 
connected, in the sense that a particle can escape from a 
trap into any line-segment (rather than just neighbour- 
ing segments as happens in reality). For intermediate 
trapping probabilities, we might expect the true distri- 
bution of first passage times to fall somewhere in be- 
tween the predictions of these two limiting cases. The 
analytical approach proceeds in a similar way for both 
approximations; we give here the details for approxima- 
tion[Tl indicating only briefly the modifications necessary 
for approximation [21 

Let us start by considering a particle which enters the 
upper lane at a random position in a tau segment of 
length N. The total sojourn time spent by the particle 
in the upper lane tu{N) is a random variable given by 
the following sum 



tu{N) = T{N) + ^ r« + ^T«(7V,). 



(28) 



Here m is the number of times the particle is reflected 
from the trap, T{N) is the time to reach the trap for 
the first time (from a random initial position), T^^'^ is 
the time spent at the trapping site on the ith visit and 
Ti^\Ni) is the time to return to the trap after escaping 
for the ith time. Under approximation 1, the particle 
is always reflected back into the same line segment and 
Ni ~ N. Averaging over stochastic histories then yields 

{tu{N))i = {T{N)) + {m + l)(Tt,) + (m)(Ti(7V)) (29) 

where we have used the independence of the various ran- 
dom variables in order to factorize the expectation values. 

Noting that the number of reflections has a geometric 
distribution given by 



P{m) = A{l-Ay 
we immediately find 



A 



(30) 



(31) 



Furthermore the time spent at the trap is exponentially 
distributed with mean 



^ A particle which enters the upper lane at a trap site is randomly 
assigned to the line segment to the left or the right. 



{Ttr) 



1 



2D + UJa 



(32) 
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Together with (|24|) and (|27p . this allows a complete de- 
termination of {tu{N))i. 

Next we note that the probability of entering the upper 
lane in a given line segment (i.e., a sequence of consecu- 
tive sites occupied by tau) is proportional to its length. 
Hence for a uniform distribution of r-particles we should 
average {tu(N))i over the distribution 



The resulting infinite sum 



(33) 



{tu)l 



(34) 



contains summands proportional to N^pt^ with a 
1,2,3 and evaluates to 



{tu)l = 



Pr 



1 



2D{l^Prf L0a{l-Pr)' 



(35) 



The only difference in approximation[2]is that the lengths 
Ni appearing in (j28p are drawn from the distribution 
P{Ni) = /9^'~^(1 — Pt)- After averaging over this distri- 
bution as well as the initial distribution one finds 



{tu)2 = 



Pt 



1 



2i:>(i 



(36) 



In similar fashion one can obtain an expression for 
{tu{N)'^) from ([28|) . The resulting calculations are 
straightforward but tedious and are outlined for conve- 
nience in the appendix. Using approximation [T] one even- 
tually obtains 



{tl)i = {iolpril + Prf + ^DuJaPAl ' pl) 
+ ADHI- prfil+'^pr+pl)] 

x{2{l-prfD^L0l} 



(37) 



Results dSlI) and ^ together with (l2T])-(l23l) finally de- 
termine the parameters of the asymptotic Gaussian form 
of P(x,i). 



3. Results for velocity and diffusion coefficient 

After the effort involved in finding the long-time be- 
haviour of P{x,t), it is comparatively straightforward 
to evaluate the integral — one simply uses Stirling's 
approximation for the factorials in the Poisson distribu- 
tion Px{yi,t) and then employs a saddle-point treatment. 
This finally yields the large deviation function of the ob- 
served average current along the lower lane 



e(yi) 



yi 

2 



yi 111 



yi 

px+ 



{pA-Pi,)x+ , Pj 



2A 



b 

2A 



(38) 



where x+ is the positive solution of the quadratic equa- 
tion 



[pA - Pb)x - yiA = 0. 



(39) 
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FIG. 8: (Color online) Comparison of analytical and simula- 
tion results for the dependence of single motor velocity (along 
microtubule) on concentration of tau. Parameter values are 
p = 0.5, D — 1.0, uja = 1.0 and cod = 0.01. Simulation results 
are for a system of size L = 10000 with a fixed realization 
of disorder for each density value; averages are taken over 
100000 time steps and 10000 histories. 



The large deviation function gives knowledge of all cu- 
mulants of Yi. In particular, we easily obtain the mean 
velocity 



Vl 



and diffusion constant 



; lim 

pPb 



Di = lim 

t^OO 



t 



^pPb+p^A. 



(40) 
(41) 

(42) 
(43) 



In fact, this is simply the expected result for the variance 
of a compound process. However, we emphasize that 
the large deviation expression (|38)) also gives complete 
information about the asymptotic behaviour of higher 
moments. 

Using (PT|) and ([^ for the bound-state occupation 
probability and its variance, we obtain explicit expres- 
sions for velocity and diffusion. Specifically, the observed 
average motor velocity along the lower lane in the pres- 
ence of a quenched distribution of r is given by (under 
approximation 1) 



p{l-prf 



i^-prr' + mpr 



(44) 



In Fig. [8] we compare the analytical predictions of both 
approximate schemes with simulation results for the 
mean velocity of a single motor. We find that both ap- 
proximations reproduce well the decrease of the average 
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FIG. 9: (Color online) Comparison of analytical and sim- 
ulation results for the dependence of single motor diffusion 
constant (along microtubule) on concentration of tau. Pa- 
rameters as in Fig. [8l 
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FIG. 10: (Color online) Crossover timescale tc = -Di/(ui)^ for 
parameter values p = 0.5, D = 1.0, uJa = 1-0 and ujd = 0.01. 
Parameters as in Fig. [H Note logarithmic scale. 



FIG. 11: (Color online) Comparison of analytical and simu- 
lation results for mean current of K-particles along the lower- 
lane (microtubule) versus concentration of tau. Motor density 
is fixed at pK ~ 0.1 and hopping rates are p = 0.5, D — 1.0, 
LOa = 1.0 and ujd ~ 0.01. Simulation results are for a system 
of size L = 1000 with a fixed realization of disorder for each 
density value; averages are taken over 10000 time steps and 
100 histories. 



B. Low motor density 

The single motor calculations of the previous subsec- 
tion can be used to give approximate results for the low- 
motor-density limit of the full model. The key observa- 
tion is that, for realistic parameters, the total unbound 
density is much smaller than the total bound density (ex- 
cept for Pr very close to unity). It is thus reasonable to 
neglect the exclusion interaction both in the upper lane 
and also for inter-lane moves. In other words we have an 
asymmetric exclusion process on the filament coupled to 
random walks (without exclusion) in the surroundings. 

In this approximation one then assumes a homoge- 
neous bound density given by 

Pb = 2pK X Pb. (45) 



velocity with increasing tau concentration and that, as 
anticipated, the simulation results lie between the two 
curves. 

Similarly the analytically-calculated diffusion constant 
Di compares well with simulation, as shown in Fig.|9l In 
particular, we see that although the system is relatively 
robust to low concentrations of tau, for high concentra- 
tions there is a dramatic increase in the variance which 
could be physiologically significant. 

In fact, it is perhaps more illuminating to consider the 
diffusion constant divided by the square of the velocity as 
shown in Fig. [TOl This defines a crossover timescale tc = 
Di/ivi)"^. For times less than diffusive fluctuations 
dominate whereas for larger times the directed motion is 
dominant. 



In other words a proportion Pf, of the total K-particlcs are 
expected to be in the lower lane (where Pb is the bound- 
state occupation probability calculated for the single mo- 
tor case). Assuming no correlations between the bound 
densities on neighbouring sites, the associated mean par- 
ticle current is given by 

(j) ^2ppKPbil-2pKPb)- (46) 

As seen in Fig. 1111 this approach is a considerable im- 
provement over the simple mean- field (i.e., neglecting all 
correlations and simply replacing Ua by the effective ab- 
sorption rate given by (U) since it implicitly takes ac- 
count of correlations between tau and unbound motor 
density. However, it still assumes an uncorrelated uni- 
form density of bound motors (only expected to be a 
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FIG. 12: (Color online) Site-dependent motor densities {bi} 
and {ui) in lower and upper lanes of the system studied in 
Fig. [TT] with pT = 0.9. The data is again averaged over 1000 
histories with a fixed realization of disorder; lines are provided 
as an aid to the eye. 



FIG. 13: (Colour online) Simulation results for variance 
aj — (j^) — (j)'^ in current of K-particles along the lower- 
lane (microtubule) versus concentration of tau. Parameters 
as in Fig. 1111 Dotted line shows expected variance for single- 
lane ASEP with density 2pKPb', peak at high tau-density is 
attributed to inter-lane fluctuations. 



good approximation when many hops along the filament 
take place between absorption/desorption events). To 
test this assumption we plot, in Fig. [TH averaged den- 
sity profiles for the rate parameters used in Fig. [TT] with 
Pr — 0.9. In fact there is significant spatial dependence 
in the bound density but the relative variation is much 
less than in the unbound density (where we find peaks 
corresponding to the trapping of particles in tau-occupied 
regions). As a quantitative comparison we remark that 
for the scenario shown in Fig. [T^] the standard deviation 
of the lower-lane densities is approximately 0.14 of their 
mean (over 1000 lattice sites) whereas the equivalent ra- 
tio for the upper lane is about 0.84. For larger p the run 
length of the molecular motors increases relative to the 
trap spacing so that the bound density becomes much 
more homogeneous while the unbound density is essen- 
tially unchanged (not shown). 

Based on the single motor results we also expect high 
concentrations of tau to significantly increase the vari- 
ance of the bound K-particle current. Indeed, this is 
observed in simulations as shown in Fig. 1131 For low 
tau-dcnsities we see a plateau where the variance is dom- 
inated by that of the ordinary single-lane exclusion pro- 
cess, which can be calculated exactly [H, [H^l- At a 
high tau concentration, a peak appears resulting from 
the inter-lane fluctuations (just as in the single motor 
case). It seems difficult to obtain the shape of this peak 
analytically since fluctuations along the lower lane and 
fluctuations between the lanes are not independent. 

V. DISCUSSION 

This work is concerned with the stochastic modelling 
of molecular-motor traffic inside cells. In particular, wc 



were motivated by experimental observations of the influ- 
ence of the linker protein tau on transport in axons [1, Q ■ 
These experiments suggest that the primary effect of tau 
proteins is to reduce significantly the binding-affinity of 
the motor proteins rather than to alter their dynamics 
on the microtubules. 

In order to model this situation, we extended an 
earlier lattice-based approach due to Klumpp and 
Lipowsky [l^. Specifically, whereas the original model 
contained only a single species of particles (representing 
kincsin motors), wc introduced a second type of parti- 
cle to represent the tau proteins. Following Klumpp and 
Lipowsky we modelled the axon in a simplified fashion 
using a one-dimensional track embedded in a cylindrical 
geometry. The dynamics of the model was defined to 
reflect the picture suggested by experiment so that the 
local presence of tau reduces the rate for kinesin particles 
to attach at a potential binding site but does not affect 
the rates for movement along the track. 

Simulation results for this model show that for a wide 
range of parameters, the presence of tau causes only a 
marginal reduction in the maximum transport along the 
axonal filament. (This finding is consistent with the very 
recent results of [5l[.) However the position of the maxi- 
mum is shifted to higher kinesin densities and for a given 
density of motor proteins one may observe a vigorous 
alteration of the flow. This result indicates that an over- 
expression of tau proteins may perturb the intracellular 
transport signiflcantly. The central aim of the present 
work was to characterize and understand this disruption 
within simple models and thus gain tentative insight into 
possible mechanisms in the real biological context. (We 
have not attempted to make quantitative predictions and 
further work would be needed to determine to what cx- 
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tent our observations are seen in real physiological con- 
ditions.) 

In this spirit, in order to analyse in more detail the 
origin of the transport disruption we further simplified 
the modelling approach by introducing a two-lane model, 
which is analytically treatable for selected, but generi- 
cal, limiting cases. For the case of the pure model, i.e., 
the case without "r-particles" , it is known that the sta- 
tionary distribution is given by a product measure [l8l |. 
However, with the introduction of r-particles, mean-field 
properties are only found in certain limits. For example, 
the case of fast diffusion of r-particles can be described 
by introducing effective binding rates for the kinesin- 
particles. 

A qualitatively different behaviour is observed if we 
consider slow diffusion of r-particles. This corresponds 
to quenched disorder on the track leading to non-trivial 
density profiles of the "K-particles" which represent the 
motors. For arbitrary model parameters this limit is dif- 
ficult to analyse. However, there are two limiting choices 
of motor rates which still yield a product measure for the 
motor distribution — these choices correspond to a space- 
dependent distribution of particles in one lane and a ho- 
mogeneous distribution in the other. In these cases one 
observes that the position of the maximal current can be 
significantly changed by the r-particles. A shift of the 
maximum to higher K-densities (for fixed r-densities) 
is also a feature of simulation results for more general 
parameters. However, when correlations cannot be ne- 
glected one also sees a decrease in the height of the max- 
imum. In these respects the simplified two-lane model 
reproduces to a large extent the properties of the full 
model. 

In the context of the biological system, it is instruc- 
tive to analyse the features of the model for a fixed den- 
sity of motor proteins and different concentrations of tau- 
proteins. Within the two-lane setup we were able to ob- 
tain analytical results for the random walk behaviour of 
a single K-particle moving among a fixed background 
of r-particles. This approach also gives approximate in- 
formation about the low-motor-density limit of the sim- 
plified model. The calculations and supporting simula- 
tions demonstrate that for a wide range of r-densities 
(with fixed motor density), the current on the filament 
remains unchanged. At high r-concentrations, however, 
one observes a drastic decrease of the system's capac- 
ity in a rather small interval of r-densities. In addition, 
the fluctuations of the current show a pronounced max- 
imum, such that the crossover-times, which characterize 
the transition from diffusive to directed motion diverge. 

It should be emphasized that these results were ob- 
tained in the context of a simplified "toy" model and cau- 
tion must be applied in drawing conclusions about real 
physiological systems. This is particularly true because 
we assumed throughout a strong reduction of the kinesin 
binding-rates, compared to the experimentally observed 
values (although this would be partly compensated by 
shorter run lengths in bidirectional transport). With 



this caveat in mind, the model results do indicate one 
possible mechanism for tau-induced transport disruption 
in neural diseases. On one hand, the robustness of the 
transport capacity for a wide range of tau-densities indi- 
cates that the system is stable against small fluctuations 
of the tau-concentration. Bearing in mind that tan pro- 
teins stabilize the network of microtubule filaments, it is 
expected that the transport capacities of nerve cells are 
only marginally influenced by the tau-concentration for 
non-pathological conditions. On the other hand, for low 
motor densities, a more extreme overexpression of tau 
can drastically reduce the transport properties of the cell 
and may trigger a cascade of events leading to a loss of 
functionality — see, e.g., [s^l for alternative mechanisms 
of such cascades. 

To conclude, we have presented a model approach to 
intracellular transport, which supports the possible im- 
portance of tau proteins for neural diseases. In future 
work the present approach should be extended, e.g., the 
structure of the nerve cell should be described in more 
detail, in order to obtain a more quantitative description 
of the relevant processes. 
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Appendix A: Mean square sojourn time in unbound 

state 

Here we summarize the intermediate steps leading to 
the result (j37]) for the mean square sojourn time in the 
upper lane. 

Starting from (j28|) one obtains (under approximation 

1) 

{U{NY),^{T{Nf) + {m.+ l){Tl) 

+ (m(m + l))(Ttr)2 + {m){Ti{Nf) 

+ (m(m-l))(Ti(7V))2 

+ 2{m+l){T{N)){T,,) 

+ 2{m){T{N)){T,{N)) 

+ 2(m(m+l))(Ti(iV))(rtr). (Al) 

In addition to quantities already introduced, we need the 
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following averages The analogous calculation under approximation 2 gives 

instead 

(-^) = ^^^^ (A2) 
(T(iV)2) ^ (A5) 

(4)2 = + Pr)' + ^DiJaPril - Pr) 

(with the latter two obtained following [1^). Finally, +4D^(1 — p^)^} 



averaging over the distribution Po(-^) of line-segment 
lengths yields 

oo 

= {ujlpr{l+ Prf +SDuJaPr{l- pI) (A7) 
+ AD''{1- Prf{l+Apr+pl)} (A8) 

x{2{l^p^fD^u:l}-\ (A9) 



x{2(l-p.ri?2^2|-i^ (AlO) 
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